PRACTICAL DESIGN OF MULTI-CHANNEL 
OVERSAMPLED WARPED COSINE-MODULATED 

FILTER BANKS 



M.I. Vashkevich*, W. Wan^^, A. A. Petrovsky* 

*Belarusian State University of Informatics and Radioelectronics, Computer Egeneering Department, 
St. P.Brovky, 6, 220013, Minsk, Belarus, email: {vashkevich, palex} @ .bsuir.by 
^Shanghai University, School of Communication and Information Engineering, 
Shanghai 200072, RR. China, email: wanwg@staff.shu.edu.cn 



Keywords: Nonuniform filter bank, optimization. 
Abstract 

A practical approach to optimal design of multichan- 
nel oversampled warped cosine-modulated filter banks 
(CMFB) is proposed. Warped CMFB is obtained by all- 
pass transformation of uniform CMFB. The paper ad- 
dresses the problems of minimization amplitude distor- 
tion and suppression of aliasing components emerged 
due to oversampling of filter bank channel signals. Pro- 
posed optimization-based design considerably reduces 
distortions of overall filter bank transfer function taking 
into account channel subsampling ratios. 

1. Introduction 

Multirate filter banks are widely used as computation- 
ally efficient and flexible building blocks for subband 
signal processing fQ. These filter banks decompose 
an input signal into its subband components and recon- 
struct the signal from the downsampled version of these 
components with little or no distortion. 
In various applications such as noise reduction, speech 
enhancement and audio coding nonuniform time-fre- 
quency representation is highly desired. A well-known 
example is approximation of critical bands of human au- 
ditory system. By using nonuniform filter banks this 
problem can be effectively solved. Another problem 
that is solved by means of nonuniform filter banks is 
an estimation of the frequency dependent reverberation 
time |6|. 

One simple way of obtaining nonuniform filter bank is 
employing allpass transform to uniform filter bank [lOj . 



Efficient nonuniform DFT polyphase filter banks were 
proposed in |[2[|5|. However, DFT-based filter banks 
produce complex- value channel signals even for real- 
value input. Therefore the subsequent subband process- 
ing becomes more sophisticated. In contrast to DFT 
filter banks allpass transformed (or warped) cosine- 
modulated filter banks (CMFB) developed in ||7| are al- 
low to avoid having complex subband signals for a real- 
value input. It is worth to mention that aliasing can- 
cellation conditions |4 | do not hold for warped CMFB 
thus all aliasing components should be suppressed by 
synthesis filter bank. 

The paper presents a practical approach to optimal de- 
sign of multichannel oversampled warped CMFB with 
low aliasing and amplitude distortions. Just like uni- 
form CMFB, analysis and synthesis filters of warped 
CMFB are obtained from one prototype filter, that re- 
sults in high design efficiency. A distinguishing charac- 
teristic of warped CMFB is that for each channel sub- 
sampling factors should be determined separately |8|. 
The practical rule of selection subsampling factors is 
also derived. 

2. Warped cosine-modulated filter bank 

2.1 Uniform CMFB: brief review 

The M-channel warped CMFB proposed in fl] based 
on uniform CMFB developed in |4|. The impulse re- 
sponses of the analysis (hk [n]) and synthesis (fk [n]) fil- 
ters are cosine-modulated versions of the prototype fil- 
ter h[n]\ 



hk[n] = 2h[n] cos 



(2/c+l)7r 
2M 



in 



1 



fk[n] = 2h[n] cos 



2M 



where < n < A^, < A: < M and i9/e = -^^ 
(-l)'f. 

The length of h[n] is assumed to be multiple of 2M, 
i.e. N = 2mM. The transfer functions of analysis and 
synthesis filters can be expressed as follows: 

(2a) 

F,{z) = alb,F{zW^i^'-') + a,blF{zW-^'^'-'^), 

(2b) 

where ak = e-^'^"^)''?, 6^ = 1^2^°"^^ ^ 
g-jTr/M jj^g superscript * denotes the complex conju- 
gation. In ^ H{z) is linear phase lowpass FIR filter 
prototype with cutoff frequency Uc ^ tt/M. 

2.2 AUpass transformed CMFB 



Allpass transformation of CMFB consists in replacing 
all the delays in uniform CMFB with causal and stable 
allpass filters: 



z 



A{z). 



(3) 



In this paper we consider first-order allpass filter 

^-1 



A{z) 



a 



1 -\- az 



a e . 



\a\ < 1, 



with frequency response A{e^^) = e^'^^^\ The phase 
response is written as 

(piuj) = —u + 2 arctan . (4) 

\ a cos uj — 1 J 

Thus, replacing all terms z~^ by by first order allpass 
filters leads to a transformation u ^{oo) of the fre- 
quency scale as shown in figure [T] 

3. Problem statement 

Figure [2] shows the block diagram of the nonuniform 
oversampled filter with overall transfer function 



Tall{z) = Tdist{z) + Talias{z), 



(5) 



where 



M-l 



Td^st{z) = Hk{z)Fk{z), 

k=0 
M-l Sk-l 

Taliasiz) =J2J2 Hk{zW'sJPk{z). 




Figure 1 : Allpass transformation of CMFB 
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Figure 2: Allpass transformation of CMFB 



Distortion transfer function Tdist{z) shows amplitude 
distortion introduced by analysis and synthesis filters, 
while aliasing distortion is described by an aliasing 
transfer function Taiias{z). In perfect reconstruction 
c^sc Taiiasie'^) = 0^nd\Tdist{e'^)\ = 1,Vcj G [0, tt]. 
Thus optimization problem for the design of the over- 
sampled warped CMFB is formulated as follows 

mm\E{u)\, VcjG[0, tt], (6) 

where the error function E{uj) is defined as 

E{uj) = \Taii{e^n?-l, uje[0,7r]. (7) 

A weighted least square objective function can be for- 
mulated to solve optimization problem ([6]), 



(8) 



u;G[-0, tt] 



k=0 1=0 



where B{uj) is nonnegative weighting function. 
It is difficult to obtain the optimal solution that mini- 
mizes ([8]). In |3j| it has been shown that an iterative 
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algorithm can be derived for the design of CMFB that 
jointly optimizes the filter prototype coefficients and 
weighted function. The design algorithm given in this 
paper adopts this approach for the minimization of ([8| to 
archive solution that is near optimal in minimax sense. 

4. The design algorithm 
4.1 Preliminaries 

It is well known that frequency response of linear phase 
FIR filter can be written as 



where 



2cos(f) 2cos(%^) ...2cos ^ — ^ 
N even 

'12 cos (cj)...2cos(i^^^) 
TV odd 



h=[h{[f\)h{[f\^l)...h{N-l)]. 

The operator [a\ denotes the largest integer smaller than 
a. The superscript ^ indicate matrix transpose. Thus the 
overall transfer function ^ can be expressed as follows 



(9) 



where 



M-1 Sk-l 
/c=0 ^=0 

and 

71 (cj, /, k) = -(l){uj + 27il/Sk) - 27r(A: + 0.5)/2M, 
72(cJ,/,fc) = -(^(cj + 27r//5'fc) + 27r(/c + 0.5)/2M. 

4.2 Optimization procedure 

The design algorithm consists of two iterative proce- 
dures: one nested within the other. The vector of coef- 
ficients h.i are iteratively optimized in inner loop proce- 
dure. In the outer loop procedure the weighted function 



B^{uj) is updated. During inner loop procedure B^{uj) 
is kept constant. Subscripts £ and /i are used to indicate 
values at the ^-th iteration of inner loop procedure and 
/i-th iteration of outer procedure, respectively. 
Using ([9]) the objective function ([8| at the ^-th iteration 
can be expressed as 

g{h,)= B^{oj)\E,{uj)f, (10) 

u;G[-0, tt] 

where 

Vr{uj) = Re(U(cj)), Vi{uj) = Im(U(a;)). 

The necessary condition for minimization of ([8) is 
equality of gradient ^ zero. Differentiating (10) 
with respect to h 

V^(h^)=2 ^ B^{uj)E^{uj)-gmd{E^{uj)), (n) 

uje[o, tt] 

where 

gTad{Ei{uj)) = 
2(hfU,(c^)h^(U,(a;) +U^(a;))h^+ 

hfU,(c^)h,(U,H +Uf (a;))h,). 

Let h£ opt = h£ + e£ denotes the optimum solution that 
sets ([TT]) to zero. Expanding in Taylor's series 

V^(h^+e^) = v^(h£) + V^^(h^)-e^ + . • • = 0. (12) 
Without taking into consideration higher order term 



in ( 12 ) the correction vector can be obtained by solv- 
ing 

= -[v'^(h,)]-i-v^(h,). (13) 
V^^ is hessian of objective function that is given by 

V'^(h^)=2 B^{uj)[gT8.d{E^{uj))x 

a;G[0, tt] (14) 

gmd{Ee{u)f ^ Ee{uj) X G(a;)], 

where 

G(c^) =2 (h^(U,(c^) + U^(c^))^(U,(c^) + Vj{uj))h^ 
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Iterative procedure starts with ho, which is a rough es- 
timation of optimum, that can be obtain using available 
FIR filter design method (windowing, etc.). At the l- 
th iteration, the correction vector obtained by solv- 



ing ( 13 ). Then coefficient vector is updated using 



-10 



inner loop procedure terminates when | |e/ 1 p < 10 
The outer loop procedure updates B^{uj) as described 
below. Let h.^^opt denote the optimum solution ob- 
tained by inner loop procedure for weighting function 
B,^{ijo). The error function E^{uo) is obtained using 



|7| and h.^^opt- Let V^{uJi) be the l-\h extremum of 
\E^{uj)\. \iv^{uji) < min[F^(a;,_i), then 
let V^^{(jJi) equal this threshold. Define the envelope 
function /3^{uj), u G [0, tt] as 



w = 



where / is chosen such that oJi < OJ < cj^+i. This enve- 
lope function shows that a larger relative weight should 
be assigned at u where \Ejj^{u) \ is large. 
The new weighting function is obtained as 



where the normalization factor is defined as 




and 6> is a positive factor that affects the convergence 
(1 < < 1.5). Iterative algorithm terminates when 

max i3^{uj)-min/3^{uj) ^ , 
max/3^(u;)+min/5/_i(u;) — ^' 

where approximate value of = 0.5 . . . 0.7. 



aliasing cancellation conditions do not hold. Therefore, 
the subsampling factors Sk in the k-th subband has to be 
chosen such that the aliasing components in each sub- 
band do not overlap with subband spectrum of input sig- 
nal. Selection of range [/l^, fuf,] should be based on 
the same principle as the mechanism of aliasing cancel- 
lation nonadjacent channels do not overlap. This means 
that the range [/l^, fuk] should include passbands of 
(A: - l)-th A:-th and (A: + l)-th subfilters. 
The general algorithm for finding the boundaries of the 
frequency band [/l^, fuk] defines as follows. The or- 
dered set 



^ ~ \ M J 



0, 1 ... M, 



(16) 



defines passband frequency of subfilters of uniform 
CMFB. Mapping ^ allows to obtain corresponding fre- 
quency edges of warped CMFB (figure [3]) 



{F^} 




{R} 



Figure 3: Transformation of frequency axis by map- 
ping Q. 

According to proposed rule the frequency band for k-th 
channel of warped CMFB expressed as 

/l. =<p(ffc"_i)/27r, /c;, = (p(F«+2)/27r, (17) 

for /c = and /c = M — 1 the following relation holds 



ho = 0, fuo = v{F^)/27r, 



(18a) 



5. Subsampling ratio selection 

In (8) the following rule is given for selection of sub- 
sampling ratios for warped CMFB 



>Sk> 



1 <nk < 



fUv, —fh 



(15) 

where fu^ fhk normalized lower and upper edge of 
the /c-th subband. 

In the uniform case /^^ and fu^ are match the passband 
of /c-th subfilter, that is possible due to mechanism of 
aliasing cancellation |j4J. However in the warped CMFB 



(18b) 

Joint use of rule (p?])-(p^ allows determining the sub- 
sampling ratios Sk for warped CMFB. 



6. Design example 

Let us consider the 2 2 -channel warped CMFB approxi- 
mating the psychoacoustic Bark scale for standard sam- 
pling frequency 16 kHz. At first it is necessary to de- 
termine warping coefficient a. According to |9| for 
/s = 16 kHz a = 0.5783. Considering rule described 
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in previous section it is possible to choose the subsam- 
pHng ratios Sk (figure |2]) such that ahasing does not af- 
fect channel signals. For instanse examine the subsam- 
pling factor S3. Necessary frequency bounds of uniform 
CMFB given bellow 



rpu 27r 

^2 — 22' 



^5 ~ 22 • 



Using ^rf\ the corresponding frequency range of 
warped CMFB can be determined as 



: 0.0122, h 



Us 



: 0.0316, 
(19) 



The subsampling factor Ss is obtained by applying ( 15 ) 



l<n3< 



fu3 



fUn -fh 



= [1.6303J = 1 ^ ns 

' >^3> 



0.0632 



0.0244 



^ [15.81J > 5'3 > [0] ^83 = 15. 

Similarly, the remaining subsampling factors for filter 
bank under consideration are selected 

Sk = {56,27,20,15,12,21,18,15, 

13, 12, 10, 9, 7, 6, 5, 5, 4, 3, 1, 1, 2, 3}. (20) 

The proposed design method was implemented using 
MATLAB on an Intel Celeron 2.8 GHz with 1 GB 
physical memory. It took nine outer loop iteration 
for the algorithm to converge (6 minutes). The fre- 
quency responses of initial and optimized filter proto- 
types are plotted in figure |4] It can be seen that pro- 
posed optimization procedure considerably minimized 
the stopband energy of filter prototype. Figure [5] shows 
the resulting magnitude frequency response of resulting 
warped CMFB. 

For a chosen subsampling ratios Sk the magnitude re- 
sponse of aliasing transfer functions for initial and opti- 
mized warped CMFB were calculated (figure [6]). It can 
be seen that aliasing distortion has the same order of 
magnitude as the stopband attenuation of the prototype 
filters in figure]?] 

The level of aliasing component which appears due 
to decimation/interpolation of channel signals can be 
shown using bifrequency system function 1 1 1 (figure [7]- 
[8]). Figure[8]reveals that optimized warped CMFB atten- 
uates aliasing component to the level of -80. . .-90 dB. 
Overall transfer functions of initial and optimized filter 
banks are given in figure [9] Overall transfer function 
of initial warped CMFB suffers from irregular distor- 
tion caused by aliasing. With optimized filter prototype 




-Optimized filter prototype 
Initial filter prototype 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

u;/27T 

Figure 4: Magnitude responses of initial (using | 3 1) and 
optimized filter prototypes of order N = 176. 




Figure 5: Magnitude response of the warped 22-channel 
CMFB (a = 0.5783). 

the ripples of overall transfer are decreased significantly 
(from 0.15 dB to 0.004 dB). Thus the design example 
shows that proposed algorithm effectively minimizes 
the overall distortion introduced by warped CMFB. 

Conclusion 

A practical method for the design of multichannel over- 
sampled warped CMFB with low level of amplitude dis- 
tortion has been proposed. Formulation of optimization 
problem and imposed constraints on overall filter bank 
transfer function are allowed to minimize amplitude dis- 
tortion. Also the rule for selection of subsampling factor 
in warped CMFB is derived. Using the proposed design 
method it is possible to obtain high quality nonuniform 
filter bank with low distortion level. 
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T^|.^^(e^°') funciton of initial warped CMFB 
T^|.^^(e''°) funciton of optimizedwarped CMFB 




Figure 6: Magnitude responses of aliasing transfer func- 
tions. 




u;/27r 



Figure 7: Magnitude bifrequency response of initial 
warped CMFB. 




UO' /2'K 



Figure 8: Magnitude bifrequency response of optimized 
warped CMFB. 
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I "'"all^®'"^ function of optimized warped CMFB 
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Figure 9: Magnitude responses of Taii{z). 
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